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We analyze the two-mode squeezed harmonic oscillator and the feth-order har- 
monic generation within the framework of Bargmann-Hilbert spaces of entire func- 
tions. For the displaced, single-mode squeezed and two- mode squeezed harmonic 
oscillators, we derive the exact, closed- form expressions for their energies and wave 
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functions. For the /cth-order harmonic generation with k > 3, our result indicates 



that it does not have eigenfunctions and is thus ill-defined in the Bargmann-Hilbert 
O" 1 space. 
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1 Introduction 



Recently there is renewed interest in formulating and solving dynamical systems involving 
harmonic modes in the framework of Bargmann-Hilbert spaces [1, 2, 3]. For example, in 
[4, 5], we applied the Bargmann-Hilbert space technique to obtain the exact solutions of 
families of quantum non-linear optical as well as spin-boson models. In [6, 7, 8, 9, 10], the 
authors applied the technique to the quantum Rabi model, a simple system describing 
the interaction of a two-level atom with a harmonic mode. A Bargmann-Hilbert space 
is a Hilbert space of entire functions introduced by Bargmann and Segal. It is a vector 
space with typical orthonormal basis A= v n = 0, 1, 2, ■ • •. Elements in the space are entire 
functions and the space is equipped with a well-defined Hermitian scalar product, 



(f,h) = j f(z)h(z)dfx(z) (1.1) 

for any two elements f(z), h(z) in the space, where d/J,(z) = -e~' 2 ' dxdy. In a Bargmann- 
Hilbert space, the harmonic creation and annihilation operators a),a can be realized as 
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a^ — >■ z, a — » j-. This realization enables one to convert the time-independent Schrodinger 
equation of a dynamical system into a differential equation. Solutions to the differential 
equation are entire functions. 

In this paper, we apply Bargmann-Hilbert spaces to analyze the two-mode squeezed 
harmonic oscillator and kth order harmonic generation. For the cases of the displaced, 
single-mode squeezed and two- mode squeezed harmonic oscillators, we derive the exact, 
closed-form expressions for their energies and wave functions. For the /cth-order harmonic 
generation with k > 3, our result shows that it does not have eigenfunctions and thus is 
ill-defined in the Bargmann-Hilbert space. 



2 Two-mode squeezed harmonic oscillator 

The Hamiltonian of the two-mode squeezed harmonic oscillator reads 

H = cj(a[ai + a\a2) + g (a\a 2 + 0102), (2.1) 

where we assume that the boson modes are degenerate with the same frequency u, and g 
is a real constant. In terms of the operators K±, K defined by 

K + = a\a\, K_ = aia 2 , K = -(a\ai + 0^2 + 1), (2.2) 

the Hamiltonian (2.1) can be written as 

H = 2u (k - i) + g {K + + K_). (2.3) 

The operators K±, K form the su(l, 1) Lie algebra. The quadratic Casimir of the algebra, 
C = K + K_ — Kq(Kq — 1), takes the particular values C — k(1 — k) in the representation 
(2.2), where k = 1/2, 1, 3/2, • • •. This is the well-known infinite-dimensional unitary 
irreducible representation of sw(l,l) known as the positive discrete series T> + (k). Thus 
the Fock-Hilbert space decomposes into the direct sum of infinite subspaces T-t K labeled 
by « = 1/2,1,3/2,- ... 

The basis state in the subspace "H K , denoted as \k,ti), n = 0, 1, 2, • • •, has the form 



A\n+2K-l( r ,\\n 

(n + 2n- l)\n\ 



and the action of K±, K in this representation is given by 



K \k, n) = {n + k)\k, n), 

K + \K,n) = yJ(n + 2K,)(n + l)\K,n + l), 

K_\k,u) = yJ(n + 2K-l)n\K,n- 1). (2.5) 



Using the Fock-Bargmann correspondence 



i d 

a 1 — > z, a — > — , 

dz 

we can show that the infinite set of monomials 



\n) 



(2.6) 



VkM) 



(n + 2n- l)\n\ 



n = 0,1,2, 



(2.7) 



form the basis in the Bargmann-Hilbert subspace associated with the representation (2.5). 
Thus the operators K±,K (2.2) have the single-variable differential realization in the 
subspace labeled by the Bargmann index k, 



d 
K = z— + K 

dz 



K 4 



K_ = z— + 2k— , K = 1/2, 1,3/2, 
dz 1 dz 



(2.8) 



By means of this differential representation (2.8), we can express the Hamiltonian (2.3) 

(i.e. (2.1)) as the 2nd-order differential operator in each Bargmann-Hilbert subspace 

labeled by k, 

( d 1\ / d 2 d\ , nn . 

Then the time-independent Schrodinger equation gives the differential equation for wave- 
function ip(z), 



d 2 d 

9 z ~f^^( z ) + 2 ( uz + 9K)-j-i>(z) 



gz + 2u [ k - - 



E 



With the substitution 



i>{z) 



e 9 



(1-A) 2 



(p(z) 



A 



-s. 



ip(z) = 0. (2.10) 



(2.11) 



where 



< 1, it follows, 



dp = l gz— + 2[uAz + gn}— + 2/twA - u - E 1 <p = 0. (2.12) 

This differential equation is exactly solvable, This is seen as follows. First of all, let 
us recall the characterization of exact solvability of a differential operator. A linear 
differential operator £ is exactly solvable if it preserves an infinite flag of finite-dimensional 
functional spaces, 

Vi C V 2 C • • • C V M C • • • , 

whose bases admit explicit analytic forms, that is there exists a sequence of finite- 
dimensional invariant subspaces Vm, M. = 1, 2, 3, • • •, 



CV M C Vm, dimVx < oo, V M = span{£i, • • • , ^dimv^}- 



In our case, we have, for any positive integer n, 

Cz n = [(2n + 2k)wA - u - E] z n + n(n + 2k- l)gz n ~ x . (2.13) 

It follows that C preserves an infinite flag of finite dimensional spaces V\ C V2 C • • • C 
Vm C ■ ■ •, with explicitly determined subspaces Vm — {1, 2^, ^ 2 , • • • , z M }, and exact solu- 
tions are polynomials in z in the Bargmann-Hilbert space. We thus seek solutions of the 
form to the differential equation (2.12), 

M 

<p(z) = \\{z- Zi ), M = 0,1,2,.-., (2.14) 

i=l 

where tp(z) = 1 for M. = 0, Ai is the degree of the polynomial and Zi are roots of the 
polynomial to be determined. Substituting into (2.12) and dividing on both sides by (p(z) 
give rise to 

M I M 9 M -1 

E + u -2kuA = gz^ Yl ^2[uAz + qk]Y^ 

i=l Z ~~ Z i j^i Z i ~ Z j i=l Z ~ Z i 

M Res 

= 2nwA + ]T ^^ £i , (2.15) 

i=l ^ z ' 1 

where Res 2=Zi are the residues of the right hand side of the first equality at the simple 

poles z — Zi, i.e. 

M 2 
Res 2=2i = gzi Y 1 - 2uAzi + 2gn. (2.16) 

^ 2j - Zj 

The left hand side (2.15) is a constant and the right hand side is a meromorphic function 
with simples poles at z — Zi. The right hand side is a constant if and only if the coefficient 
of all the residues at the simple poles are vanishing. We thus obtain the energies 



E = -u + 



■2M • 2 ( k - - 



uA, (2.17) 



and the system of algebraic equations satisfied by the roots Zi, 

5]^— + -A + - = 0, i = l,2,---,M. (2.18) 

jjki Z i Z j 9 z i 

The corresponding wavefunctions are given by 

M 



-2(i-A)* 

i=i 



4,(z) = e—V- A)z l[(z-z i ). (2.19) 



3 /cth-order harmonic generation 

The Hamiltonian of the A;th-order harmonic generation reads 



,U k J. n k 



H = outfa + g (a T ) fe + a k , (3.1) 



where k = 1,2, ••• is any positive integer, and g is a real constant. The k = 1 and 
k = 2 cases of (3.1) give the Hamiltonians of the displaced and single- mode squeezed 
harmonic oscillators, respectively These two oscillator models can be solved by the single- 
mode Bogoliubov transformation [11]. For k > 3, (3.1) gives models with higher order 
harmonic generation. 

Introduce the operators Q±, Qq in terms of the harmonic mode, 

Q +=w (a,)i - Q -=m ak - Qo= K a ' a+ D- (X2) 

Then in terms of Q±,Qo, the Hamiltonian (3.1) can be written as 

H=ku [q q - 1) + gVk* (Q + + Q_) . (3.3) 

It can be shown [4] that the operators Q±, Qq form a polynomial algebra of degree k — 1, 
defined by the commutation relations 



where 



[Qo,Q±] = ±Q±, 




[Q + ,Q-] = (fc) (Q o )-0 W (Qo-l), 




w ^=-n(^ + i-i) + nC**- 


i 



(3.4) 



(3.5) 



is a A; </l -order polynomial in Q . The Casimir operator of the algebra is given by 

C = Q_Q + + 0«(Q O ) = Q+Q- + (fe) (Qo - 1). (3.6) 

For k = 1 and k = 2, (3.4) reduces to the Heisenberg and sw(l, 1) algebras, respec- 
tively. Thus, the algebra (3.4) can be viewed as polynomial deformation of su(l, 1) and 
Heisenberg Lie algebras. 

The realization (3.2) provides a unitary irreducible representation of the polynomial 
algebra, which for k = 2 reduces to the well-known positive discrete series of su(l, 1). In 
the realization, the Casimir (3.6) takes the particular value, 

°-a(^-p)- (37) 

If we use q to label the basis states of this representation in the Fock-Hilbert space Tib, 
then it can be shown that q takes k values, 

1 jfc + l 2k + I (jfc-l)jfc + l 

For k = 2, then C = ^ and g equals to \,\, as expected. Thus the boson realization (3.2) 
corresponds to the infinite dimensional unitary representation with particular q values 



1 (fc-l)fc+l 

(3.8), and the Fock space Tib decomposes into the direct sum 7i b = 7i b ® • • • ®l-i b k 

i (fe-i)fc+i 

of k irreducible components % b , ■ ■■,'H b k 

The basis state \q,n), n — 0, 1, • • • , in the irreducible representation space T-L b is then 
given by [4] 

a' 



,1k(n+q—k) 



\q,n) = 

^j[k[n + q-±- 2 

The action Qo,Q± in this representation reads 

Qo\Q,n) = (q + n)\q,n), 

* / ik-1 

Q+\Q,n) = [[[n + q + 



= |0), n = 0,l,2, 



(3.9) 



f,2 



\q,n+ 1), 



„ i \ A/ (i-l)A; + l\ , 
Q-\q,n) = 1.1 n + ? o 1 l^™ -1 )- 



j=i 



A- 2 



(3.10) 



By using the Fock-Bargmann correspondence (2.6), we can make the following associ- 
ation 

\q, n) — ► *,,„(*) = = = . (3.11) 



,A;(n + g-^)j! 

It can then be shown that in each subspace "H 9 of the Bargmann-Hilbert space with 
basis vectors ty qtn (z), the operators Q±, Qo (3.2) are realized by single- variable differential 
operators 



Qo 

Q + 

Q- 



d 

z~r + q, 

az 



{\fk) k ' 

2 - W n(4+, a- 1 ** 1 ' 



i=i 



ob 



A- 2 



(3.12) 



We remark that there is no singularity in the differential operator expression for Q-. This 
is because the z~ l term disappears in the expansion of the product in Q_ thanks to the 

fact that rijLi (q ~ k 2 + ) — ^ ^ or an ^ ne a ll° we< i Q values. 

Using the differential realization (3.12) we can equivalently write (3.3) (i.e. (3.1)) 
as the /cth-order single-variable differential operator in each Bargmann-Hilbert subspace 
labelled by index q, 



H = kuJ ( z i +q -h) +9 



z + k* z 



k -1 



k ( d 



(j-i)fc + r 

A; 2 



(3.13) 



Thus the time-independent Schrodinger equation for the model yields 



z k z- 1 



k ( d 



(j-i)fc + r 

A; 2 



gz + ku lz— +q-^j-E},=(). « 3.1 I) 



This is a k-th order differential equation of Fuchs' type. Solutions to this equation 
must be analytic in the whole complex plane if E belongs to the spectrum of H . In other 
words, we are seeking solution of the form 

oo 

il>(z) = '£K n (E)z n , (3.15) 

n=0 

which converges in the entire complex plane, i.e. solution ip(z) which is entire. 
Substituting (3.15) into (3.14), we obtain the 3-step recurrence relation, 

K 1 (E) + A K (E) = 0, 

K n+1 (E) + A n K n (E) + B n K n ^(E) = 0, n > 1, (3.16) 



where 



A r . 



^(n + q-^-i;) 



g^nLf^ + i^-^^^' 



i=i V'" ' ' H fc 2 



Bn ~ 1 k nk I , 1 , (j-l)k+l\- ( 3 - 17 ) 

« Ilj=i [n + l + q- y—j£ — J 

The coefficients A n , B n have the behavior when n — > oo, 

A n ~an a , B n ~bn p (3.18) 

with 

a=-^, a = -k + l, 6=-!, f3 = -k. (3.19) 

Thus the asymptotic structure of solutions to the n > 1 part of (3.16) depends on the 
Newton-Puiseux diagram formed with the points Po(0, 0), Pi(l, — k + 1), P2(2, — /c) [12]. 
Let a be the slope of Po-P an d r the slope of P1P2 so that a = a and r = (3 — a. 
Applying the Perron-Kreuser theorem (i.e. Theorem 2.3 of [12]), we have 

k = 1 case: the displaced harmonic oscillator, a = 0, r = —1, that is o > t. In this 
case, the truly 3-term part (i.e. the n > 1 part) of the recurrence relation (3.16) has two 
linearly independent solutions K n) i,K Hy 2 for which, when n — > 00 

-fr'n.i gk k ~ v K Hy2 ku 

So X™ n = K n2 is a minimal solution of the truly 3-term part of (3.16) with k — 1. The 
corresponding infinite power series solution is generated by substituting K™ m for the K^s 
in (3.15) and converges in the whole complex plane, i.e. it is entire. 

k = 2 case: the squeezed harmonic oscillator, a = — l,r = —1, and so a = t = a. 
The characteristic equation of the n > 1 part of (3.16), t 2 = 0, has two equal solutions 
ti — ^2 — 0. Then all solutions of (3.16) behave similarly as n — > 00, viz, 

lim sup (\K n \n\)™ =0 (3.21) 



for all non-trivial solutions of the 2nd equation of (3.16) with k = 2. The zero limit 
means that the n > 1 part (i.e. the truly 3-term part) of the recurrence (3.16) possesses 
a minimal solution K™ m and the corresponding infinite power series expansion, obtained 
by substituting K™ m for the K n 's in (3.15), converges in the whole complex plane, i.e. it 
is entire. 



k > 3 case : a = —k + 1, r = —1, and thus point Pi lies below the line segment PqP\ in 
the Newton-Puiseux diagram. Then 



i_ 



Jim sup (|tf n | («!)*)"=-= (3.22) 



for all non-trivial solutions of the 2nd equation of (3.16). This indicates that solutions to 
the truly 3-term part of the recurrence (3.16) with k > 3 are dominant, and the corre- 
sponding infinite series expansion (3.15) has a finite radius of convergence proportional to 
vP. It follows that there does not exist solution to (3.14) with k > 3 which is entire. We 
thus conclude that the /cth-order harmonic generation model with k > 3 does not have 
eigenfunctions (and is ill-defined) in the Bargmann-Hilbert space. 

By the Pincherle theorem, i.e. Theorem 1.1 of [12], the ratios of successive elements 
of the minimal solution sequences K™ m for the k = 1, 2 cases are expressible in terms of 
infinite continued fractions. Proceeding in the direction of increasing n, we find 

(3.23) 



T/'mm 

itn — ~ — " 

T/'mm 

n 


B n +\ 


B n +2 -Dn+3 

A n +2~ A n+ %— 


which for n = gives 






T/'min 

R - l 


A,- 


B2 B 3 


' T/'min 
n 


A 2 - A 3 - 


ry'min 

Note that the ratio R = J- min involves KV} m , although the 



(3.24) 



expression is obtained from the truly 3-term part of (3.16), i.e the recurrence (3.16) for 
n > 1. However, for single-ended sequences such as those appearing in the infinite power 

jy-min 

series expansion (3.15), the ratio R$ = * jw of the first two terms of a minimal solution 
is unambiguously fixed by the n = part (i.e. the first equation) of the recurrence (3.16), 
namely, 



w (? - h - £) 



Ro ~ ~ Ao ~ -gv-inUb + i-^' (3 ' 25) 

In general, the Rq computed from (3.24) is not the same as that from (3.25) (i.e. (3.24) 
and (3.25) are not both satisfied) for arbitrary values of recurrence coefficients A n and 
B n . As a result, general solutions to the recurrence (3.16) are dominant and are usually 
generated by simple forward recursion from a given value of Kq. Physical meaningful 
solutions are those that are entire in the Bargmann-Hilbert spaces. They can be obtained 
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if E can be adjusted so that equations (3.24) and (3.25) are both satisfied. Then the 
resulting solution sequence K n (E) will be purely minimal and the power series expansion 
(3.15) will converge in the whole complex plane. 

Therefore, if we define the function F(E) = Rq + Aq with Rq given by the continued 
fraction in (3.24), then the zeros of F(E) correspond to the points in the parameter space 
where the condition (3.25) is satisfied. In other words, F(E) = yields the eigenvalue 
equation, which may be solved for E by standard nonlinear root-search techniques. Only 
for the denumerable infinite values of E which are the roots of F(E) = 0, do we get entire 
solutions of the differential equations. 

As a matter of fact, the spectra for the k — 1,2 cases can be determined explicitly. 
As will be seen in the next two subsections, the infinite power series in (3.15) actu- 
ally truncates for these two cases, so that their solutions are given by polynomials in 
Bargmann-Hilbert spaces. 

3.1 Displaced harmonic oscillator 

The displaced harmonic oscillator is the k — 1 special case of the /cth-order harmonic 
generation. By (3.14), the time-independent Schrodinger equation in the Bargmann- 
Hilbert space of analytic functions reads 

(cjz + g)^ + (gz-E)^ = 0. (3.26) 

(3.27) 



With the substitution 




rj}{ z ) = e - 9z/ "<f)(z 


the above differential equation reduces to 




.<<" + »>£ -(*+£). 



[z) = 0. (3.28) 

This differential equation is exactly solvable, and exact solutions are polynomial of the 
form 

A/" 

(/>(z) = H(z-Zi), AT = 0,1,2,..-, (3.29) 

where <fi(z) = 1 for M = 0, M is the degree of the polynomial and z\ are the roots of the 
polynomial to be determined. Following the procedure similar to that in the last section, 
we obtain the energies of the system, 

E = uUr-£\, (3.30) 

and the set of algebraic equations determining the roots Z{, uzi + g = 0, i = 1, 2, • • • , Af. 
It follows that Z{ = —- and the solution (3.29) has the form 

9\ ( . 9\ M 



%=\ 



n * + S = * + S ■ (**: 



Thus the wavefunction of the model is given by 



N 



9\ 

u> 



(3.32) 



These expressions for the energies and wavefunction agree with those in [1] by a different 
approach. 



3.2 Single-mode squeezed harmonic oscillator 

The Hamiltonian of the single-mode squeezed harmonic oscillator is given by 

H = coa'a + g 



f^) 2 + a 2 



(3.33) 



which corresponds to the k = 2 case of the /cth-order harmonic generation Hamiltonian 
(3.1). In the Bargmann-Hilbert space, the time-independent Schrodinger equation reads 



d 2 

4gz—ij(z) + (2uz + 



d 
dz 



i/>(z) + 



<]■-- 2u;[q--) -E 



ij(z) = 0, 



(3.34) 



where q — |, | are the Bargmann index of su(l,l). This equation is the k = 2 special 
case of the fcth order differential equation (3.14). 
With the substitution 



where 



2.9 



tp(z) = e 4 s^ ' z (p(z), 
< 1, it follows, 



n 



i 



CO 2 '' 



(3.35) 



d 2 



d 



Agz-— + YluVLz + 8gq] — + 2quVt u - E } <p = 0. 

dz 2 dz 2 



(3.36) 



This differential equation is exactly solvable, and exact solutions are polynomials in z 
which automatically are entire functions in the Bargmann-Hilbert space. We thus seek 
solution of the form, 



ip(z) 



M 

no*-*) 



M = 0,l,2, 



(3.37) 



where tp(z) = 1 for M. = 0, M. is the degree of the polynomial solution and z% are the 
roots of the polynomial to be determined. Following the procedure similar to that in the 
last section, we obtain the energy eigenvalues, 



E 



1 



2,M + 2( g -i 



UjQ,, 



and the set of algebraic equations which determine the roots Zi 

2 



M 

E 



&i Zi z i 



to „ 2o 
2g Zi 



1,2,---,M. 



(3.38) 



(3.39) 
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The corresponding wavefunctions are 

M 



^(z) = e--9 (l ~ il)z \{{z-z i ). (3.40) 
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